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ABSTRACT 

MAXI J1836—194 is a Galactic black hole candidate X-ray binary that was discovered 
in 2011 when it went into outburst. In this paper, we present the full radio monitoring 
of this system during its ‘failed’ outburst, in which the source did not complete a full 
set of state changes, only transitioning as far as the hard intermediate state. Obser¬ 
vations with the Karl G. Jansky Very Large Array (VLA) and Australia Telescope 
Compact Array (ATCA) show that the jet properties changed significantly during 
the outburst. The VLA observations detected linearly polarised emission at a level of 
~1% early in the outburst, increasing to ~3% as the outburst peaked. High-resolution 
images with the Very Long Baseline Array (VLBA) show a ~15mas jet along the po¬ 
sition angle —21 =b 2°, in agreement with the electric vector position angle found from 
our polarisation results (—21 ±4°), implying that the magnetic field is perpendicular 
to the jet. Astrometric observations suggest that the system required an asymmetric 
natal kick to explain its observed space velocity. Comparing quasi-simultaneous X-ray 
monitoring with the 5 GHz VLA observations from the 2011 outburst shows an un¬ 
usually steep hard-state radio/X-ray correlation of Lr oc L]^ 8±0 ' 2 , where Lr and Lx 
denote the radio and X-ray luminosities, respectively. With ATCA and Swift monitor¬ 
ing of the source during a period of re-brightening in 2012, we show that the system 
lay on the same steep correlation. Due to the low inclination of this system, we then 
investigate the possibility that the observed correlation may have been steepened by 
variable Doppler boosting. 

Key words: X-rays: binaries - radio continuum: stars - stars: individual (MAXI 
J1836-194) - ISM: jets and outflows 
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1 INTRODUCTION 


Relativistic jets are observed from actively-accreting black 
holes on all scales, from stellar-mass black holes in low- 
mass X-ray binary (LMXB) systems to supermassive black 
holes at the centres of galaxies. There is a clear but poorly 
understood relationship between the morphology and spec- 
trum of the jets and the s tructure of the accretion flow (e.g. 
iFenderl l20ld ; [Gallol l2010| j . LMXBs provide ideal laborato¬ 
ries to study jets and their connection to the accretion flow 
as they evolve through their full duty cycles on humanly- 
observable timescales. Spending the majority of their life¬ 
times in a low-luminosity quiescent state, LMXBs occasion¬ 
ally go through a much brighter outburst phase (the re¬ 
sult of increased accretion onto the central compact object) 
where the jet properties change significantly as the system 
transitions through a number of canonical X-ray states (see 
IBelloni et al.ll201ll for a full review of X-ray spectral states). 

In a typical outburst, LMXBs are initially observed in 
the hard state (HS), which is associated with persistent ra¬ 
dio emission from a steady, partially self-absorbed, compact 
jet | Dhawan et al.l 120001 : ICorbel et al.l l2000l : IStirling et al.l 
12001 ) that has a flat or inverted , optically-th ick radio spec¬ 
trum (a > 0, where S v oc //*; IFenderl l200ll h As the out¬ 
burst progresses, the accretion rate increases and the sys¬ 
tem may transit through a number of intermediate states 
as it approaches the soft state. During the transition to 
the full soft state the jet emission transforms into discrete, 
bright, relativis tically-moving knots with an optically -thin 
radio spectrum (ICorbel et al]|2004l : lFender et al.ll2004l h The 
comp act jet is quenched by at least 2 . 5 orders of magni¬ 
tude (iFender et all Il999al : ICoriat et all l201ll : iRussell et all 
12011 3) as the jet spectral break moves to lower frequen¬ 
cies iCorbel et aHl2013bl: IRussell et al.ll2013bl ; I Russell et all 


l2014al l and particle accelerati on seemingly occurs at differ¬ 
ent locations within the jet (IRussell et al.l l2014a| j . At the 
end of the outburst, the X-ray luminosity decreases and the 
source transitions back to the hard state. In this reverse 
state transition the compact jet is gradually re-established 
(first switching on in the radio band and then the in¬ 
frared / optical band; M iller-Jones et all 120121 : ICorbel et all 


l2013bl : [Kalemci et al 


2013l l. indicating the potential build 


up of particle accelerating structures increasingly closer to 
the black hole. 

The processes that are responsible for launching the 
jet outflow are poorly understood. However, there exists a 
well-studied relationship between the radio and X-ray lumi¬ 
nosities in hard state LMXBs, which ha s been observed for 
both individual sources (e.g. GX 339—4; lHannikainen et all 
1 19981 : ICorbel et all 1200,1 2013b|)_and the entire sample of 
BH X-ray binaries ie.g. IGallo et al.|l2003l) . Although recent 
work has cast some do ubt over the universality of the rela¬ 
tion (IGallo et al .1120 1 -fl ) . it is generally described by two dis¬ 
tinct power-law tracks in the radio/X-ray luminosity plane, 
an upper track and a lower, steep e r, ‘radio-quiet’ trac k 
(IGallo et all 120121 : ICoriat et all l201ll : ICorbel et all l2013bh . 
The track followed by a given source does not appear to 
be connected to any property of the object (such as spin; 
ISoleri fe Fendeill201ll 'l and there is evidence for some indi¬ 
vidual sources transitioni ng between the two t racks as their 
X-rav luminosity fade s (ICoriat et all 1201 ll : Ijonker et all 
I 2 OI 2 I : iRatti et alj|20ljl . It has also been suggested that this 


bimodal relationship may also exist in AGN, supporting the 
theory that the mechanisms responsible for extracting en- 
ergy from the acc retion flow are similar across the mass scale 
iKing et al.ll2013h . 

High resolution, very long baseline interferometry 
(VLBI) observations can reveal important information 
about the changing properties of the jet during an outburst. 
Directly imaging the jet can determine the orientation of the 
jet axis, the bulk motion and structure of the jet, and how 
these change with the evolving accretion flow. While discrete 
ejecta have been observed in a number of systems in out- 


burst (e.e. Mirabel &; Rodrieuez 1994J; 

Tineav et al.lll995l; 

Hiellming & Ruoenl 119951: Fender et al. 

1999q), the hard 


state compact jet has on ly been spatially re solved in two 
syste ms, GRS 1915+ 105 dDhawan et al.ll200d l and Cygnus 
X—1 (IStirling et al.|[200ll ~). both of which are not standard 
LMXBs. High-cadence observations of the discrete ejecta 
can determine their proper motions. Extrapolating back¬ 
wards in time, the moment of launching may then be coupled 
to changes within the accretion f low to identify the mech - 
anisms responsible for launching dMiller-Jones et ai1l2012l) . 
Astrometric observations of a system over time may also al¬ 
low for the proper motion of the system to be calculated, 
which can then be used to constrain the formation mecha- 
nism of black holes (for full discussion, see lMiller-Jonedl2014l . 
and references therein). 


1.1 MAXI J1836-194 


MAXI J1836—194 was discovered at the onset of its outburst 
in August 2011, triggering an intensive multiwavelength ob¬ 
serving campaign. Classified as a black h ole candidate due to 
its radio and X-ray prop e rties dstrohmaver fc Smithl 2011 


ay properties (btrolimayer & bmitn 
Miljer-Jones_etj^l201lj: IRussell et all 2011al : Reis et ail 
20121 : IRussell et alj |2014bh . MAXI J1836—194 transitioned 


from the hard state (spectral states are defined b y X-ray 
spectral and timing properties: iFerrigno et akll2012 ) to the 
hard intermediate state (HIMS) on 2011 September 11. 
The outburst ‘failed’ t o enter the full soft sta te, remain¬ 
ing in the HIMS (see iBrocksopp et all 120041 for further 
discussion on failed outbursts), reaching its softest point 
on September 16, following which it transitioned back to 
the hard state on September 28 and faded towards qui¬ 
escence. In 2012 M arch, the system underwent a period 
of renewed activity ferimm et all 120121 : lYang et al .1 1201 (3 : 
iGrebene v et al.l 201 ,'jl) befo re returning to quiescence by 2012 
July 05 ( Yang et al .11201 21)1). This low-inc lination system (be¬ 
tween 4° and 15°: IRussell et a~L i]2014bD is between 4 and 
10 kpc away (IRussell et all 2014al l and has proven itself 
to be an ideal system in which to study th e evolution of 
the radio jet due to its bri ght radio emission (IRussell et all 
l2013bl : IRussell et al.|[2014al 'l . In this paper we present results 
from an intensive VLA, ATCA and VLBA radio monitor¬ 
ing campaign of MAXI J1836—194, which began during the 
early phase of the outburst (2011 September 3) and con¬ 
tinued until the source became sun constrained (our last 
radio observation occurred on 2011 December 3). We also 
present ATCA and quasi-simultaneous Swift monitoring of 
the source during its 2012 re-flare. In particular, we will dis¬ 
cuss the evolution of the polarised radio jet, the jet morphol¬ 
ogy, and the VLBI astrometry during the 2011 outburst, as 
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well as the hard state radio/X-ray correlation as the system 
faded towards quiescence. 

2 OBSERVATIONS 

2.1 Radio observations 

2.1.1 VLA 

From 2011 September 3 through December 3, we conducted 
an intensive monitoring campaign on MAXI J1836—194 with 
the Karl G. Jansky Very Large Array (VLA), under project 
code VLA/10B-218. We have previously presented the six 
epochs of VLA observations with simultaneous milli metre, 
mid-/near-infrared, optical and X-ray data llRussell et al.l 
l2013bl : iRussell et al.l I2014al) . Here, we report the full 18 
epochs of VLA radio observations covering the rise, tran¬ 
sition to and from the HIMS, and subsequent decay of the 
2011 outburst of MAXI J1836—194. The VLA observations 
were taken every 2-5 days during the brighter phase of 
the outburst (September 3 - September 30) and every 1 
- 2 weeks as the outburst decayed, until the final obser¬ 
vation on 2011 December 3. The array was in its most- 
extended A configuration until September 12, in the A—>D 
move-configuration until September 18, and thereafter in the 
most-compact D configuration. 

The observations were taken over a wide range of fre¬ 
quencies (between 1 and 43 GHz) with an integration time of 

3 seconds for all observations prior to November 11, follow¬ 
ing which the data were taken with an integration time of 
5 seconds. The lowest-frequency 1-2 GHz observations were 
recorded in two 512-MHz basebands, each containing eight 
64-MHz sub-bands made up of sixty-four 1-MHz channels. 
All other frequencies were recorded in two 1024-MHz base¬ 
bands comprised of eight 128-MHz sub-bands made up of 
sixty-four 2-MHz channels. Following standard data reduc¬ 
tion procedures with in the Common Astr onomy Software 
Application (CASA; iMcMullin et al.1 120071) the data were 
corrected for shadowing, instrumental issues, and radio fre¬ 
quency interference. 3C286 was used for bandpass and am¬ 
plitude calibration, J1407+2827 was used as a polarisation 
leakage calibrator, and phase calibration was performed us¬ 
ing the nearby (~ 7° away) compact source J1820—2528. 
Calibrated data were then imaged on a per-baseband basis 
using the multi-frequency synthesis algorithm within CASA 
and subjected to multiple rounds of self-calibration. To de¬ 
termine the flux density of the source, we fitted a point 
source to the target in the image plane (Stokes 7) and stan¬ 
dard VLA systematic errors were added. Stokes Q and U 
were also measured at the position of the peak source flux in 
the lower-frequency (<10GHz) data (Table[T]). At no epoch 
was the source observed to be significantly extended. 

A nearby source is located to the south-west of 
MAXI J1836—194 (^(Uaway; see Figure[lJ inset). This un¬ 
related, steep-spectrum source was not aligned with the jet 
axis and was not variable over the duration of the out¬ 
burst. At times when the VLA was in its most compact 
D-configuration (for all observations after 2011 September 
18), our low-frequency (<4GHz) VLA observations could 
not resolve the target from the confusing source. Therefore, 
where unresolved, our observations of MAXI J1836—194 are 
contaminated by the nearby source. We must determine the 



Frequency (Hz) 


Figure 1 . VLA spectrum of MAXI J1836—194 (blue points) 
and the nearby unrelated source (red points) from 2011 Novem¬ 
ber 1. At frequencies where the nearby object could not be re¬ 
solved from our target (<4GHz) we subtracted the expected 
MAXI J1836—194 brightness from the combined flux densities 
to approximate the emission from the unrelated source. Inset: 
November 1 VLA 5 GHz image of MAXI J1836—194 (black cross), 
showing the confusing source to the south west of the target. Con¬ 
tour levels are at ±(y/2) n times the rms noise (0.1 mjy), where 
n = —3,3,4,5... (dashed contours represent negative values). 
At low frequencies (<4GHz) when the VLA was in the com¬ 
pact D configuration, MAXI ,11836—194 was unresolved from the 
nearby unrelated object. Fitting S„ against v with a power law, 
we then estimate the spectral index of the confusing source to be 
—1.81±0.02. 


spectrum of the unrelated object to estimate its contribu¬ 
tion and then subtract from necessary observations. Using 
an epoch when MAXI J1836—194 was faint (2011 Novem¬ 
ber 01), we measured the brightness of the confusing source 
(Figure [TJ. At frequencies where we were unable to re¬ 
solve the two sources from each other, we subtracted the 
expected MAXI J1836—194 brightness (determined by fit¬ 
ting a power law to the higher frequencies) from the total 
spectrum to give an approximation of the unrelated source. 
Fitting the measured radio flux densities, Sj,, against fre¬ 
quency, v, with a power law, we then estimate the spectrum 
of the nearby source and MAXI J1836—194. The spectrum 
of MAXI J1836—194 was well fit by a single power law at all 
epochs. 


2.1.2 VLB A 

MAXI J1836—194 was observed five times with the VLBA 
during the outburst, for between 4 and 6 hours per epoch, 
under project code BM339 (Table [5]). We observed in dual 
polarisation mode, with 64 MHz of bandwidth per polarisa¬ 
tion. In the first four epochs, we made use of the dual 13/4- 
cm recording mode to split the bandwidth equally between 
two frequency bands centred at 8.4 and 2.25 GHz. However, 
since our VLA monitoring indicated the source flux density 
to be < 1 mjy at the time of the final epoch, we observed 
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Table 1. Sample VLA flux densities of MAXI J1836—194. a is the radio spectral index, la errors are uncertainties on the fitted source 
parameters. Stokes Q and U are given before systematic errors are added. The full table is available online. 


Date 

(UT) 

MJD 

Frequency 

(GHz) 

Flux density 
(mjy) 

Q 

(mJy beam -1 ) 

u 

(mjy beam -1 ) 

a 

2012 Sep 03 

55807.1 

4.60 

27.2T0.3 

0.24T0.05 

-0.29i0.06 

0.64i0.03 



7.90 

38.4T0.5 

0.67T0.06 

-0.49i0.07 


2011 Sep 05 

55809.1 

5.00 

23.2T0.3 

0.08T0.02 

-0.08i0.02 

0.84i0.04 



7.45 

32.5T0.4 

0.18T0.02 

-0.16i0.02 


2011 Sep 07 

55811.2 

1.50 

18.1T0.4 

- 

- 

0.51i0.01 



5.00 

28.6T0.7 

0.25T0.07 

-0.13i0.07 




7.45 

39.3T0.9 

0.58T0.09 

-0.08i0.08 




20.80 

68±2 

- 

- 




25.90 

79±2 

- 

- 




32.02 

88±5 

- 

- 



Table 2. Elliptical Gaussian fits (in the uu-plane) to the 2.3 and 8.4 GHz VLBA images of MAXI J1836—194, showing the flux density, 
beam size, major axis, minor axis and the position angle for each VLBA epoch. The radio flux brightened during the X-ray peak of the 
outburst and the major axis size remained steady at 8.4 GHz during our observations. 


Date 

2011 

(UT) 

MJD 

Freq. 

(GHz) 

Flux 

density 

(mjy) 

Beam 

size 

(mas) 

Major 

axis 

(/xas) 

Minor 

axis 

(//as) 

Position 

Angle 

(°) 

Sep 04 

55808.12 i 0.06 

8.42 

30.8 i 0.2 

2.5 x 0.8 

584 i 52 

188 i 27 

-17.0 i 4.6 



2.27 

15.9 i 0.3 

16.5 x 4.8 

2950 i 200 

2180 i 1380 

90.0 i 0.1 

Sep 09 

55813.08 i 0.08 

8.42 

38.5 i 0.2 

2.6 x 0.8 

699 i 33 

140 i 45 

-25.7 i 3.3 



2.27 

24.6 i 0.3 

12.2 x 3.4 

6810 i 320 

2180 i 170 

-17.5 i 2.5 

Sep 19 

55823.06 i 0.08 

8.42 

29.3 i 0.1 

2.4 x 0.8 

685 i 32 

101 i 46 

-25.6 i 3.0 



2.27 

28.9 i 0.2 

14.3 x 4.4 

4870 i 310 

2830 i 110 

-18.9 i 5.7 

Oct 12 

55846.02 i 0.10 

8.42 

5.0 i 0.1 

2.3 x 0.8 

833 i 123 

Oi 395 

-11.8 i 4.5 



2.27 

3.4 i 0.2 

19.6 x 0.8 

9850 i 1750 

2 i 2730 

-15.5 i 5.5 

Nov 18 

55883.88 i 0.10 

8.42 

0.36 i 0.09 

2.8 x 1.0 

- 

- 

- 


at 8.4 GHz only on November 18, using the full 64 MHz of 
available bandwidth. 

The observations were phase-referenced, using the clos¬ 
est available extragalactic calibrator source, J 1832-2039, 
from the VLBA Calibrator Survey (VCS; iBeaslev et all 
l2002ll . whose assumed position (from the most up-to-date 
global astrometric solution available at the time of the 
observation^]) was 18 h 32 m 11 s .0465, -20°39'48".202 (1.6° 
from the target source). The phase referencing cycle time 
was 3 minutes, spending 2 minutes on the target source and 
1 minute on the phase reference calibrator in each cycle. 
Every seventh observation of the target was substituted 
for a 1-minute observation of an astrometric check source. 
We used two different check sources, alternating between 
J1845-2200, taken from the Fourth VLBA Calibrator Sur¬ 
vey (VCS-4: iPetrov et al .1 [200(3 ). and J1825 -1718, from the 
Second VLBA Calibrator Survey (VCS-2; iFomalont et al.l 
I2003T) . 30 minutes at the start and end of each observation 
were used for geodetic blocks, observing several bright extra- 
galactic calibrators over a range of elevations. This improved 
the accuracy of the phase referencing process by allowing us 
to solve for unmodelled tropospheric delays and clock errors. 

The data were correlated using the VLBA-DiFX 
implementat i on of the software correlator developed by 
iDeller et al.l (|2007h . and reduced using the Astronomical 

1 http://astrogeo.org/vlbi/solutions/rfc_2011c/ 


Image Processing System (aips; iGreise 3 120031 '), according 
to standard procedures. Data with elevations < 15° were 
flagged, to prevent the increased atmospheric phase fluctu¬ 
ations at low elevations from affecting the astrometric ac¬ 
curacy. Scatter-broadening affected the data at 2.25 GHz, 
reducing the amplitude measured on the longer baselines, 
and we were unable to calibrate the station at Mauna Kea 
(MK). Phase solutions for MAXI J1836-194 were derived 
by fringe-fitting the phase reference calibrator, using as a 
model the best available image of the calibrator following 
self-calibration. 

For the first four epochs, the data were imaged with 
robust weighting as a compromise between sensitivity and 
resolution. For the first three epochs, phase-only self¬ 
calibration was performed down to a solution interval of 
5 minutes. These phase corrections were applied, before a 
single round of amplitude and phase calibration was used 
to adjust the overall amplitude gains of the individual an¬ 
tennas. By the fourth epoch, the source was significantly 
fainter, and we were only able to perform a single round of 
phase-only self calibration on a timescale of 1 hour. For the 
faint, fifth epoch, we imaged the data with natural weight¬ 
ing to maximise the sensitivity, but the source was too faint 
for self-calibration. 

Since the final images showed the source to be 
marginally resolved in the first f our epochs, we exported the 
self-calibrated data to Difmap llShepherdlll997l j. and fit the 
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source with a single elliptical Gaussian in the uv-plane. Since 
Difmap does not provide uncertainties on fitted model pa¬ 
rameters, we used the best-fitting values as initial estimates 
for the AIPS task uvfit, which was used to determine the 
final fitted source sizes. 

2.1.3 ATCA 

The Australia Telescope Compact Array (ATCA) observed 
MAXI J1836—194 four times during its 2011 outburst and 
five times in 2012 (Table [3j. All observations were taken 
at 5.5 and 9.0 GHz with a bandwidth of 2 GHz at each fre¬ 
quency. The ATCA data were redu ced using standar d rou¬ 
tines in ATNF’s MIRIAD package rtSault etal.ll 19951 ). Pri¬ 
mary flux calibration for all observations was done using 
1934—638. Phase calibration was carried out with 1817—254 
during the 2011 September 05, 20, and 21 observations. 
1908—201 was used for all other observations, except on 2012 
April 08 where 1829—207 was used. Due to a combination of 
antenna shadowing in the compact H75 configuration and a 
missing antenna, only one baseline was usable for the 2011 
September 20 and 21 observations. Therefore, for these two 
epochs a point source was fit to the calibrated data in the uv- 
plane. For all other observations, calibrated data were then 
exported to Difmap for imaging and flux densities were ob¬ 
tained by fitting a point source model to the target in the 
image plane (Table 0. Due to the compact array configura¬ 
tion on 2011 September 20 and 21, our ATCA observations 
were unable to resolve MAXI J1836—194 from a nearby un¬ 
related object. The flux density of the nearby source was 
subtracted from our measured flux density to estimate the 
flux density of MAXI J1836—194 in our ATCA observations 
(see Section [2.1.11 and Figure [T] for further discussion). 

The radio counterpart of MAXI J1836—194 was de¬ 
tected for all 2011 observations at a position coincident with 
those of the X-ray and optical counterparts. During our 2012 
observations, the source was detected only on April 08 (dur¬ 
ing the minor re-flare). When detected, the source was found 
to be unresolved within the ATCA beam. 

2.2 X-ray observations 

2.2.1 RXTE 

The outburst of MAXI J1836—194 was monitored with the 
Proportional Counter Array (PC A) on board the Rossi X - 
ray Timing Explorer (RXTE- see ljahoda et alii 19961 . 1200614 . 
A total of 74 pointed observations sample this outburst, 
for a period of about 92 days. The RXTE data were re¬ 
duced using the Heasoft software package v6.8, following 
the standard steps described in the RXTE cookbook. To 
estimate the source inten s ity a nd a spectral colour, we 
followed lAltamirano et akl (120081 ): we used the 16 second 
time-resolution Standard 2 mode of all active Proportional 
Counter Units (PCUs, all layers) in each observation. The 
hard colour was defined as the count rate ratio (6-16 
keV) /(2-6 keV) and the intensity was defined as total count 
rate in the 2.0-20.0 keV band (Figure [3]). We subtracted 
the background and performed deadtime corrections on a 
per-PCU basis. To take into account response differences 
between PCUs, we normalised the intensity and hard colour 
by the Crab before estimating the average intensity and 


Table 3. ATCA flux densities of MAXI J1836—194. Due to the 
presence of a nearby unrelated source (within the ATCA beam, 
Figure [l] inset), for the observations in which the array was in the 
compact H75 configuration the interpolated flux density of the un¬ 
related source (determined from our VLA observations; Figure[TJ 
was subtracted to estimate the flux density of MAXI J1836—194 
(see Section l2.1.1ll . At other epochs, the more extended array con¬ 
figuration meant the two sources were resolved from each other. 
Upper limits are 3<r. 


Date 

(UT) 

MJD 

Array 

config. 

Freq. 

(GHz) 

Flux 

density 

(mjy) 

2011 Sep 05 

55809.45 

6B 

5.5 

28±2 




9.0 

36±4 

2011 Sep 20 

55824.22 

H75 

5.5 

45±4 a 




9.0 

50±3 a 

2011 Sep 21 

55825.21 

H75 

5.5 

39±4 a 




9.0 

44±3 a 

2011 Dec 01 

55897.27 

1.5D 

5.5 

0.4T0.1 




9.0 

0.24T0.05 

2012 Jan 13 

55939.92 

6A 

5.5 

^0.09 




9.0 

0.27T0.07 

2012 Feb 03 

55960.05 

6A 

5.5 

^0.11 




9.0 

^0.21 

2012 Apr 08 

56024.94 

H168 

5.5 

0.64T0.05 




9.0 

0.68T0.04 

2012 Jun 24 

56102.73 

6D 

5.5 

^0.23 




9.0 

^0.28 

2012 Aug 24 

56163.30 

6A 

5.5 

^0.07 




9.0 

^0.09 


Preliminary results for the 2012 April 08 epoch were reported 
by ^joriat_et ii al i (2012). 

a The flux density of the nearby source was subtracted from our 
measured flux density (see Section 12. 1.1 1 . 

hard colour on a per-observation basis. Due to the low col- 
umn density along the line-of-sight (riH = 3.0 x 10 21 cm" 2 ; 
[Russell et al.ll2014al ) and the spectrum being dominated by 
a power-law component during the hard state decay, the 
3-9 keV band during the decay phase of the outburst was 
unaffected by absorption. 

Due to Sun constraints, no simultaneous RXTE data 
were available to coincide with the 2011 December 3 VLA 
data (the final RXTE epoch was 2011 November 25). We 
estimate the X-ray intensity on December 3 by fitting the 
observed exponential decay of the source following its tran¬ 
sition back to the hard state (see Figure O, and extrapo¬ 
lating to the required date (error bars were conservatively 
estimated from the uncertainty of the exponential fit). 

2.3 Swift XRT 

The Swift X-ray telescope (XRT) monitored 
MAXI J1836—194 during its period of renewed activ¬ 
ity in 2012. From the HEASARC public archives, we 
retrieved XRT observations taken on 2012 April 06 and 
09, June 15 and 29 (which were closest in time to our 
April 08 and June 24 ATCA observations). The 2012 April 
observations (characterised by a higher count rate) were 
taken in windowed-timing mode, while those from 2012 
June (when the source was fainter) were in photon-counting 
mode. We extracted the spectral files and their associated 
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Figure 2. Radio and X-ray lightcurves of MAXI J1836—194 during the VLA, VLB A, ATCA and RXTE monitoring. The grey shaded 
region represents the HIMS (also denoted by HIMS) and HS is the hard state, as defined bv lFerrigno et all d2Q12T) . The top panel shows 
the flux densities of the VLA, VLBA and ATCA monitoring. The second panel gives the RXTE 2-20 keV lightcurve (normalised to the 
Crab) and the third panel shows the evolution of the radio spectral index, a, where S u oc i/ a (calculated from the ATCA and VLA 
observations, see Section Em . The spectral index flattens during the HIMS, and then becomes more inverted during the decay. 


background, response and ancillary res ponse files with th e 
online XRT data product generator (lEvans et al.l l2009h . 
The X-ray spectra were then fit with t he X-R ay Sp ectral 
Fitting Package (XSPEC) version 12.7 dArnaudlll996l l. We 
modelled the spectra in the 0.4-10 keV band with a simple 
disk-blackbody plus power-law, absorbed b y a (fixed) 
line-of -sight Galactic riH = 3.0 x 10 21 cm ” 2 (IRussell et al.l 
l2014al l and free intrinsic absorption. More complex fits with 
physical Comptonisation models or irradiated disk models 
did not improve the fit due to the low signal-to-noise ratio 
of our data. Moreover, our primary interest was in the 


flux from the power-law component in the 3-9 keV band, 
rather than the details of the soft component. Therefore, 
a disk-blackbody plus power-law model was adequate. On 
2012 April 06 and 09, a faint disk-blackbody component was 
significantly detected with a temperature kT m — 0 . 16^0 02 
keV (from the combined spectra), but the 3-9 keV flux 
was entirely dominated by the power-law component 
(r = 1.65l°;J°) and was unaffected by absorption (total 
hh < 4 x 10 21 cm -2 ). On 2012 June 15, the spectrum was 
best fit by a simple power law (T = 1.98lg'g2). On 2012 
June 29, the source was no longer significantly detected. We 
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log(Radio flux density (mJy)) 

Figure 3. Hardness-intensity diagram of MAXI J 1836-194. The 
RXTE X-ray hardness (6-16 keV/2-6 keV) and intensity (2-20 
keV) are normalised to the Crab. The 5 GHz VLA and 5.5 GHz 
ATCA radio observations are shown as filled circles, with the 
marker size and colour representing the radio flux density. The 
white triangles represent the timing of the VLBI observations 
(includ ing the two CVN and EVN observations from lYang et all 
l2012al. see Section l.'i.U . VI through V5 refer to the VLBA obser¬ 
vational epochs, Cl is the CVN observation and El is the EVN 
observation. Our high-cadence radio observations show the full 
evolution of the compact jet during the failed outburst. 


therefore estimated a 90% confidence limit to its net count 
rate from the number of counts detected in the source and 
back ground regions (using Bayesian statistics; iKraft et al.l 
1 199lli . We converted this count rate upper limit to a flux 
upper limit assuming the same spectral model that we fit 
to the 2012 June 15 spectrum. 


3 RESULTS 

3.1 Source evolution 

With our VLA, VLBA, and ATCA monitoring, we ob¬ 
served the brightening and decay of the compact jet of 
MAXI J1836—194 during its 2011 outburst (Figure [5] top 
panel). The VLA, VLBA and ATCA radio observations (Ta¬ 
bles mm and [3l respectively) show that the flux density of 
the system increased during the bright, early hard state and 
HIMS, before the outburst failed. The VLA radio spectrum 
(which was initially inverted in the hard state) flattened as 
the system softened and transitioned to the HIMS (Figure^ 
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Figure 4. Polarisation parameters of MAXI J1836—194 during 
the early part of the outburst, before the VLA shifted to a more 
compact configuration and the target could not be resolved from 
the nearby confusing source. The unshaded region represents the 
hard state (HS) and the shaded represents the HIMS. Top panel: 
Stokes I at 5 and 7.5 GHz (red and blue points respectively). Sec¬ 
ond panel: linear polarisation (P) at 5 and 7.5 GHz (red and blue 
points respectively), where P= yjQ 2 + U 2 . Third panel: frac¬ 
tional polarisation (FP) at 5 and 7.5 GHz (red and blue points 
respectively), where FP= 100 x P/I. Fourth panel: the rotation 
measure ( RM ), the black dashed line is a weighted mean of the 
RM values. Fifth panel: the intrinsic electric vector position an¬ 
gle ( EVPA ; black circles), derived from fitting the observed 5 and 
7.5GHz polarisation angles on a per-subband basis against A 2 . 
The red and blue points represent the 5 and 7.5 GHz polarisation 
angles (per-baseband), respectively. We see low levels of fractional 
polarisation from MAXI J1836—194 during the early hard state 
and HIMS, a steady RM (of |RM| ^ 32 rad m —2 ) and weighted 
mean EVPA of —21 d= 4° (fifth panel, black dashed line). 
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bottom panel). Never entering the soft state, the source tran¬ 
sitioned back to the hard state on September 28 and the 
radio emission faded by ~2 orders of magnitude and the 
spectrum became more inverted during the decay towards 
quiescence. 

The X-ray evolution of a black hole X-ray binary is often 
visualised with a hardness intensity diagram (HID), which 
shows the variation of X-ray intensity with X-ray spectral 
shape (where harder spectra are dominated by the X-ray 
power-law component a nd the softer emiss ion by the black- 
body disk emission; e.g. lFender et al]|2004l . Figure [3] shows 
how the flux density of the VLA and ATCA observations 
varied with X-ray intensity and spectral shape, and indi¬ 
cates the timing of the VLA, ATCA and VLBA radio ob¬ 
servations. 


3.2 Polarisation of the radio jet 

Linear polarisation was significantly detected (> 3<r) for the 
first eight VLA epochs at both 5 and 7.5 GHz (see Figured 
second panel). The linear polarisation ( P ), fractional po¬ 
larisation (FP) and polarisation position angle (PA) were 
derived from the measured flux densities of the Stokes Q 
and U images, where P= Q 2 + U 2 , FP= 100 x P/I and 
PA= 0.5 arctan (U/Q). Following the September 18 epoch 
the target faded and the VLA shifted to a more compact 
configuration and, due to the presence of the polarised con¬ 
fusing source, we were not able to place any constraints on 
the source polarisation. 

Faraday rotation is the rotation of the plane of polari¬ 
sation as a function of wavelength (A), due to birefringence 
of the interstellar and local medium. The rotation measure, 
RM, is given by RM oc n e B^dl (where n e is the thermal 
electron density along the line of sight, B m is the magnetic 
field along the line of sight, and d is the source distance). 
The intrinsic electric vector position angle ( EVPA ) of the 
source is related to the wavelength and polarisation posi¬ 
tion angle by EVPA=PA—RM\ 2 . Therefore, the RM and 
EVPA can be determined from the 5 and 7.5 GHz observa¬ 
tions by linearly fitting the PA for all 16 subbands (ensuring 
we do not suffer from an nn ambiguity). Our results show 
that the RM and EVPA remained constant within errors 
throughout the outburst and we find that \RM\ ^ 32 rad 
m -2 and EVPA= —21 ±4° (Figure[T] fourth and fifth panel, 
respectively). 

3.3 VLBA radio morphology 

Our VLA, ATCA and VLBA radio observations show similar 
flux densities of MAXI J1836—194 during its outburst. Also, 
contemporaneous VLBI observations of MAXI J1836—194 
with the CVN (on 2011 O ctober 10) and the EVN (on 
2011 October 17) reported bv lYang et al.l d2012al l. show sim¬ 
ilar flux densities to our contempraneous radio observations, 
suggesting that the high resolution observations do not re¬ 
solve out any significant extended emission. 

The compact jet was not directly resolved in the individ¬ 
ual images, except on September 19 at 2.3 GHz. Model fit¬ 
ting in the wu-plane shows the source to be slightly extended 
at all frequencies, and in all but the last epoch (Table [2|, as 
confirmed during the deconvolution process while imaging. 


As a check on the quality of the fits, the fitted Gaussian was 
subtracted from the uw-data, and the data were re-imaged. 
The resulting residual images were noise-like, with the ex¬ 
ception of the 2.3 GHz September 19 observation, for which 
some structure remained at the level of 1.0 mjy beam -1 , sug¬ 
gesting that a simple elliptical Gaussian was not a good rep¬ 
resentation and that the source structure was more complex 
(Figure [SJ). However, since it can account for >96% of the 
emission at that epoch, we take the fitted parameters as a 
good approximation to the source structure in this epoch, 
together with a faint residual jet component separated by 
~15mas from the core, along the position angle —29 ± 3°. 
Excluding the September 04 2.27 GHz observation, which 
was taken early in the outburst with relatively low resolu¬ 
tion, model fitting in the uv-plane gave a consistent position 
angle, with a weighted mean of —21 ±2°. Hereafter, we take 
the weighted mean as indicative of the observed jet axis. 
The fitted major axis of the source (which was smaller at 
higher-freq uency, corresponding to smaller scales within the 
jet; see also lRushton et aDl2012 1 did not change significantly 
during the outburst (Table [2]). 


3.4 Astrometry 

To constrain the proper motion of the system, we measured 
the source position in each of the 8.4 GHz images prior to 
performing any self-calibration (Tabled. In all cases, we fit¬ 
ted the source in the image plane with an elliptical Gaussian, 
for consistency with the known source structure. We esti¬ 
mated systematics from the scatter on the two check source 
positions over the five epochs, weighted by their relative dis¬ 
tance to the phase calibrator. Over the two months of our 
observing campaign, the fitted source position changed by 
0.76±0.09mas in right ascension (R..A.) and 0.9±0.2mas in 
declination (Dec.). To increase the size o f our data se t , we in - 
cluded the 8.3-GHz position derived bv lYang et alj (|2012al f 
from their Chinese VLBI Network (CVN) observations. Al¬ 
though the CVN observations used the same phase reference 
source, J1832-2039, they used a less recent global astromet¬ 
ric solution to determine the assumed position of that cal¬ 
ibrator, so we corrected for the shift in calibrator position 
before adding the CVN point to our astrometric sample. 

We also attempt to correct for the motion of the opti¬ 
cal depth r = 1 surface along the jet axis, by calculating 
the expected core shift as a fraction of the comp onent of 
proper motion along the VLBI jet axis. Hein z (I2006T ) present 
a model for the expected core shift in a compact jet. Assum¬ 
ing equipartition between the particles and the magnetic 
field and following their assumptions, we calculate the ex¬ 
pected core shift at each epoch. We find a cumulative shift 
that is <15% of the total positional shift along the jet axis, 
which we include (on a point by point basis) in the uncer¬ 
tainties of our proper motion fit. 

Using jackknife resampling with a linear function to fit 
the motion in both co-ordinates, we derived proper motions 
in R.A. and Dec. of 

Ua cos 5 = — 2.3 ± 0.6 mas year -1 (1) 

Us — —6.1 i 1.0 mas year -1 , (2) 
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Figure 5. VLBA images of MAXI J1836-194. Top panels are at 8.4 GHz, and bottom panels are at 2.3 GHz. Contour levels are 
dz2 n \/2 times the rms noise level indicated in the top left corner of each image, where n = —2, 2,3,4,.... The colour scale shows the 
image flux density in units of mJybeam -1 , except for the final 8.4-GHz image (2011 November 18), for which the units are ^Jy beam -1 . 
The white crosses in the upper panels mark the predicted core position from the proper motion fit performed in Section 13.41 The source 
is unresolved except in the 2.3-GHz image of 2011 September 19, where a ~15mas extension is seen along the position angle —29 d= 3°. 


Table 4. VLBA positions of MAXI J1836—194 at 8.4 GHz at each epoch. We see the fitted position change by 0.764=0.09 mas in right 
ascension (R.A.) and 0.9=b0.2mas in declination (Dec.). 


Date 

2011 

(UT) 

MJD 

Freq. 

(GHz) 

R.A. 

Dec. 

Sep 04 

55808.12 ± 0.06 

8.42 

18:35:43.4445732(8) 

-19:19:10.48441(3) 

Sep 09 

55813.08 ± 0.08 

8.42 

18:35:43.4445671(5) 

-19:19:10.48438(2) 

Sep 19 

55823.06 ± 0.08 

8.42 

18:35:43.4445587(4) 

-19:19:10.48470(2) 

Oct 12 

55846.02 ± 0.10 

8.42 

18:35:43.444564(1) 

-19:19:10.48520(4) 

Nov 18 

55883.88 ± 0.10 

8.42 

18:35:43.444519(6) 

-19:19:10.4853(2) 


with the reference position 

R.A. = 18 h 35 m 43!444575 ± 0.000003 (3) 

Dec. = —19°19 m 10'.'4843 ± 0.0002 (4) 

at the reference date of MJD 55800.0 (Tableland Figure|H|). 


4 DISCUSSION 

4.1 The evolving compact jet 

During the initial radio observations of MAXI J1836—194, 
the system was in the hard X-ray spectral state. In this 
state, the radio spectrum was very inverted (Figure [5] bot¬ 
tom panel) and linear polarisation was observed at lev¬ 
els of a few percent (Figure [4]), significantly lower than 
the maximum expected theoretical levels for a perfectly 
ordered magnetic field. Synchrotron theory predicts linear 
polarisation at levels of up to ~10% for optically-thick 
synchrotron emission and up to ~72% for optically-thin 


synchrotro n emission in the presence of an ordered mag¬ 
netic field llPa cholczvk fc Swihartll 19671: IScheu er fc Williamsl 


1 19681 : [Rvbicki fc Lightman 19791 : L ongairl 1201 il l. H owever, 
these levels have rarely been observed in LMXBs (iFenderl 
i200df . For example, like MAXI .11836—194, linear polarisa¬ 
tion at levels of a few percent was obser ved in the hard 
states of V404 Cyg llHan fc Hiellmin jpl 9921 ) and GX 339—4 
llCorbel et, al.1 l200fll ). The low polarisation implies either 
that the m agnetic field was intrins ically disordered on large 
scales (e.g. IBrocksopp et alj|2007t l. possibly due to multi- 
ple polarised comp onents (which are polarised differently; 
IStirli n g et al.l 120041). a rotation of the polarisation angle 
within the jet ( Bland ford _fc_K diiiul 1979) or Faraday rota¬ 
tion within the jet (see Brocksopp et al.1 120131 for further 
discussion). 


As the system transitioned to the HIMS, the radio emis¬ 
sion brightened, the radio spectrum flattened, and the linear 
polarisation remained at a few percent (although it did in¬ 
crease marginally). During this time, the polarisation spec¬ 
trum steepened slightly (Figure U second panel), which the 
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Figure 6. The measured positions of MAXI J 1836—194 as a func¬ 
tion of time (positions given in Table[2]l. Here, the motion in R.A. 
and Dec. are shown in the upper and lower panels, respectively. 
VLBA measure ments are shown in red, and the 8.3-GHz CVN 
measurement of lYang et al.l 120 1 2al ) is shown in blue. We include 
the systematic uncertainty determined from the scatter in the 
check source positions, as well as the uncertainty due to the core 
shift of the jet in our error bars (applied on a point by point ba¬ 
sis). These high-resolution observations show the proper motion 
of the system. 


evolution of FP shows is not just due to the change in the 
spectral index of Stokes /. While our observations do not 
allow for a clear interpretation of this result, we speculate 
that this steepening may have been due to the magnetic 
field becoming intrinsically more ordered further along the 
jet (away from the compact object), or less ordered closer to 
the compact object. The outburst then failed in the HIMS 
and the radio spectrum became more inverted as the source 
transitioned back to the hard state. 

During the initial hard state observations of 
MAXI J1836—194, the radio spectrum was very inverted, 
with a spectral index of a ~ 0.6 — 0.8 (determined from the 
Stokes I). While X-ray binaries in the hard state typically ex- 


(e.g. Fender 200ll: McClintock et al. 

2009): Russell et al.1 

2013al; BrocksoDD et al.l 20131: van der Horst et al.l 2013!: 

Curran et al.| 2014|), inversions of a ~ 

■> 0.7 were observed 


during the decay phase of th e 2010-2011 outburst of 
GX 339—4 (ICorbel et al.l l2013al ). The radio spectrum of 
MAXI J1836—194 then flattened (to a ~ 0.1 — 0.2) as the 
system transitioned into the HIMS. Our results show that 
during the transition phases of the outburst, the radio 
spectrum appeared to track the X-ray hardness (Figure Q. 
However, during the hard state decay the spectrum of the 
radio jet continued to evolve while the X-ray hardness 
remained steady. During this time, the jet spectral break 
was observed to change by an order of magnitude indicating 
that the jet proper ti es are still c h angin g significantly 
( [Russell et al] l2013bl ; iKalemci et al.l 1 20131 : iRussell et al.l 
l2014al) . 

A flattening of the radio spectrum as the source softens 
has been observed in some (but not all) other systems. The 
shape of the spectrum is thought to be dependent on the 
geometry of the compact jet and the flow of material along 
the jet (e.g. the bulk flow velocity alon g the jet and shocks 
and turbulence within the compact je t; Blandford & Konigl 


and turbulence within the compact jet; Jblandrord & tvonigl 

ll979l ; [Falcke fc BiermaniJl995l ; lKaiseil2006l ; |Pe , er & Casella 
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Figure 7. The VLA spectral index (a) and RXTE /PCA hard¬ 
ness ratio (6-16 keV/2-6keV) of MAXI J1836—194 during its 2011 
outburst, a (left hand axis) is represented by the black points 
and X-ray hardness (right axis) is shown as the red squares. The 
grey shaded region represents the HIMS. During the HIMS and 
transition back to the HS a appeared to track the X-ray hard¬ 
ness. Following the transition back to the hard state the X-ray 
hardness remained steady while the jet continued to evolve. 


12004 1 Jamil' et al. I l2010l : iMalzaclIioTl 120141) . At the begin¬ 
ning of its 2003 outburst, the radio spectrum of H 1743—322 
showed a similar evolution to MAXI J1836—194, where the 
inverted spectrum of the compact jet flattened as the sys- 
tem softened during its transition to the intermediate state 
llMcClintock et al]|2009h . MAXI J1836—194 reached its soft¬ 
est point in the HIMS, following which the jet spectrum be¬ 
came more inverted a s the system transit ioned back to the 
hard state. GX 339-4 (ICorbel et al.|[2013al ) showed a similar 
evolution where the radio spectrum became more inverted 
following the transition back to the hard state (after the 
compact jet was re-established). However, a radio spectrum 
that steepen s as the source softens is not observed in a ll 
LMXBs (e.g. Ivan der Horst et al . 2013i : ICurran et ai.ll2l)l I i. 
From this small sample, the difference in the evolution of the 
jet spectrum between systems does not appear to be related 
to any physical parameter of the system such as inclination 
angle and black hole mass. While it appears that the evolu¬ 
tion of the radio spectrum must be driven by some property 
of the accretion flow (that may also be driving X-ray hard¬ 
ness), the surrounding environment, or the magnetic field 
(that may affect the jet flow and shape), our observations 
do not allow us to identify what processes may be driving 
the evolution. 

Interestingly, while the flux density of the target 
changed, the fitted size of the compact jet at 8 G Hz re¬ 
maine d relatively constant (Table 0 • According to iHeind 
(1200(4) . the distance of the radio core along the jet axis, zo, 
is expected to scale with the observed flux density, where 
zo oc S^ 17 . The size of the compact jet, d core , is expected 
to scale linearly with the d istance of the radio core from the 
origin dParagi et al.ll2013l ). Therefore, the fitted size of the 
compact jet should relate to the flux density of the system, 
where d co re oc sV 17 ■ In our VLBA observations we there¬ 
fore expect to see the jet change size by a factor of ~2.5 at 
8 GHz. However, our results show that the jet size remained 
steady during the outburst (even as the jet faded signifi¬ 
cantly), which could be due to physical properties of the jet 
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changing as it faded (such as changes to the opening angle, 
jet confinement region etc.). 

The VLA observations show that the RM remained 
approximately constant throughout the observations (Fig¬ 
ure U), where \RM\ <32 radm -2 (to lcr). Thi s result 
is consistent with values found by I Taylor et al.l (l2009l l 
for radio sources along a similar line of sight (where 
— 100 < RM < 100 radm -2 ). Faraday rotation of the lin¬ 
early polarised synchrotron emission is related to the ther¬ 
mal electron density and magnetic field along the line of 
sight. Therefore, the stable RM suggests that the rotation 
is not dominated by local material (to the source). The low 
value of the RM suggests that Faraday rotation is most likely 
not responsible for the low levels of observed linear polari¬ 
sation. 


4.2 Jet axis alignment 

We find an agreement between the calculated EVPA 
(—21 ± 4°; Section 13.21) and the observed VLBI jet axis 
(—21 ±2°; Section [33]). The total jet spectrum is composed 
of contributions from both the optically-thick and optically- 
thin regi ons of each of the indivi d ual synchrotron com¬ 
pone nts (iBlandford fe Konigll 1 19791 : iHiellming fe Johnstonl 
11988]) . Fractional polarisation is considerably higher in the 
optically-thin region of each synchrotron component than 
the optically-thick region; hence, it is expected that the 
observed polarisation is dominated by contributions from 
the optically-th i n reg ion of each synchrotron component 
dZdziarski et alj 12014 ) and therefore that the EVPA will 
be aligned perpendicular to the magnetic field. This im¬ 
plies that the magnetic field is perpendicular to the jet 
axis in MAXI J1836—194. A perpendicular magnetic field 
may arise fr om shock compression of the magnetic field 
Cl.ajn gl HI), knots or hotspots in the jet dDreher et alj 
|1987|), or a helical magnetic field st r uctur e dominating the 
observed emission (M arscher et alj 12008]) . a consequence 
of the rotation of the accretion disk from which the jets 


may be launched (Koide et alj 2002 

Komissarov et al. 20071; 

iMarscher et al.1 20o4" Gomez et al. 

2008lb A similar align- 


ment between jet axis and EVPA has also b een observed 
in th e black hole X-r ay binaries GX 339—4 (ICorbel et alj 
l200d) GRO J1655-40 (fll annikai nen et al.ll200Q]) and Cygnus 
X—1 dRussell fe Shahbadl20lJ . although this measurement 
is model-dependent). Polarisation studies of BL Lac sources 
have shown a bimodal distribution of the EVPA relative to 
the jet axis, where the majority of sources sh ow an align¬ 
ment between the EVPA and the jet axis (e.g. I Cohen et alj 
120141 and references therein). If an alignment between the 
EVPA and jet axis is common to black hole X-ray binaries, 
high resolution images would not be required to determine 
the orientation of the jets. 


4.3 Proper motion 

From this single outburst we calculate the proper motion 
of the system by linearly fitting its VLBI positions in 
both co-ordinates over time (see Section and Figure [6]), 
which can then be used to constrain the formation mecha¬ 
nism of the black hole. Black holes in LMXB systems are 
thought to form in two ways (|Frver &; Kalogerall20Qll h ei¬ 
ther a massive star collapsing directly into a black hole (e.g. 



Figure 8. Peculiar velocity (u pec ) of MAXI J1836—194 as a func¬ 
tion of source distance (in kpc). The black solid line is u pe c, the 
dotted lines are the error bars from uncertainties in proper motion 
and systemic velocity, and the unshade d regions denote the s ource 
distance limit of between 4 and 10 kpc dRussell et all2014al) . The 
area above the horizontal dashed line represents the requirement 
of an asymmetric supernova kick. Our results suggest that an 
asymmetric supernova kick is required to explain the peculiar ve¬ 
locity of the system. 


iMirabel &; Rodriguedl2003l ; lDhawan et ai]|2007lh or delayed 
formation in a supernova, as fallback on to the neutron star 
of materi al ejected during th e explosion creates the black 
hole (e.g. IMirabel et al.ll200l] ). Following a supernova, the 
centre of mass of the ejected material continues moving with 
the velocity of the progenitor immediately prior to the explo¬ 
sion. The centre of mass of the binary system then recoils 
in the opposite direction a nd is constrai ned to lie in the 
orbital plane (Blaauw kick: lBlaauwlll96l] ). In the presence 
of asymmetries in the supernova explosion, a further asym¬ 
metric kick, which need not be i n the orbital plane, may 


1995|; Porteeies Zwart & Yungelson 1998|; |Lai 

et al. 2001: 

Gualandris et al.1 120051: IWillcms et al.1 20051: 

Fragos et al.1 

2009]). For further discussion on formation see 

Miller-Jones 


and references therein. 


Wi th our calculated pro per motion, source distance (4- 
10 kpc; iRussell et al.l 12 014 all . and s ystemic radial velocity 


(61±15 km s 1 ; Russell et al.ll2014b] ). we can determine the 


peculiar velocity of this system. We calculate the heliocen- 
tr ic Galactic space velocity com ponents using the formalism 
of I Johnson fe Soderbloml dl987 f assuming the standard so¬ 
lar motion measured by lSchonrich et all d2010b . We assume 
circular rotation about the Galactic Centre in the plane of 
the disk, a flat rotation curve with a tangential velocity of 
240 km s - 1 , and a Galacto centric distance for the Sun of 
8.34 kpc (jReid et alj|2014bh . Defining the peculiar velocity 
as the difference between the measured three-dimensional 
space velocity and the expect ed motion of the sy stem due 
to Galactic rotation (see, e.g. iDhawan et 311200?] ). we find 
that the peculiar velocity is > 70kms -1 (Figure [8]). Since 
stellar velocity dispersion in the disk gives a typical pecu- 
l iar velocity of < 45kms -1 even for the oldest M-type stars 
llMignardll2000f ). and typical Blaauw kicks in LMXBs (while 
system parameter depe ndent) have a maximu m recoil veloc¬ 
ity of ~ 70 kms -1 (e.g. iNelemans et al.lll999f) . an asymmet¬ 
ric natal kick was most likely required in MAXI J1836—194. 

This is just the sixth system in which constraints have 


© 2015 RAS, MNRAS 000.111 [Til 










































































































12 T. D. Russell et al. 


7 10 21 

N 

Ph 

7 

03 

bp 

QJ 

X io 20 

N 

DP 

O 

LO 

'g 10 19 
£ 

I 

.2 

ii 

* 10 18 


Figure 9. The hard-state radio/X-ray correlation of 
MAXI J1836—194. The white triang les and squares are rep¬ 
resentative data from GX 339—4 (ICorbel et al] l2013bl) and 
H1743—322 (jCoriat et al.ll201lh . respectively. MAXI J1836—194 
data from the 2011 outburst decay are shown by the filled green 
and red circles for the allowable distance limits of 4 and 10 
kpc, respectively, while the open red and green circles show the 
monitoring from its minor re-brightening in 2012. The observed 
radio/X-ray correlation for this system is significantly steeper 
than either of the standard tracks shown by the representative 
data. 
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bee n placed on the for mation mechanism of the black hole 
(see iMi ller-Jones I2014I . and references therein). While mea¬ 
suring the distribution of black hole kick velocities can con¬ 
strain the formation mechanism of black holes, these obser¬ 
vations must be carried out in a quiescent or hard state, 
which are faint (high resolution observations of LMXBs are 
generally taken during full outburst). Therefore, with so few 
systems having accurate constraints, there is relatively little 
observational data on black hole kicks and observations such 
as these are required to make further progress. 


4.4 Radio/X-ray correlation 


The decay phase of the outburst allows us to probe the hard- 
state radio/X-ray correlation of MAXI J1836—194. In Fig- 
ure[9]we plot the observed hard state radio/X-ray correlation 
of MAXI J1836—194 wit h representative dat a from the up¬ 
per branch (GX 339 —4; Corbel e t al.ll2013bl ) and the lower 
branch (H1743— 322: ICoriat et al.ll201ll l. For comparison, we 
calculate th e luminosities of MAX I J18 36—194 over t h e sam e 
energies as ICorbel et al.l (l2013blJ and ICoriat et al.l d201 ill . 
During the decay, we find a steeper than usual radio/X- 
ray correlation, with an index of 1.8 ± 0.2. This correlation 
index does not appear to be dependent on frequency (e.g. 
taking any band of our radio data produces a similar index) 
and therefore cannot be explained by the changing radio 
spectral in dex. Comparing this result with expected rela¬ 
tions from Plotkin et al. (I2012i j also shows that the steep 
correlation cannot be explained by models where the X- 
ray emission originates either as inverse Compton from a 
corona, or is dominated by optically thin synchrotron radi¬ 
ation from the jet. We note that our 2011 X-ray observa¬ 
tions only span ^1 order of magnitude in X-ray luminos¬ 
ity (L x ) and deviations from the expected relationship have 
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Figure 10. Variable boosting of Lr. Top panel: the hard-state 
radio/X-ray correlation of MAXI J1836—194 (red points), with 
the de-boosted, modelled intrinsic data (blue points). Lower 
panel: V required to shift L|<. intrinsic to Lr observed i as a function 
of X-ray lumin osity assuming a dis tance of 10 kpc and an incli¬ 
nation of 10° llRussell et al.lEoi4bh . Here, log(r) oc k/ x , where 
lx = Lx/LEdd an< l K ~ 100—120 (depending on the normalisa¬ 
tion). Here, the solid blue and dashed blue lines show results from 
the minimum and maximum normalisations used to account for 
the observed scatter about the lower track. The results show that 
variable Doppler boosting of the radio emission could account for 
the unusually steep radio/X-ray correlation in this system. 


been found when the X-ray luminosity range is roughly one 
order of magnitude (thought to be due to changes in the 
spectral inde x, electron cooling, additional X-ray emitting 
regions, etc; ICorbel et al.l l2013bl) . However, if we include 
ATCA and Swift monitoring of the source during (and fol¬ 
lowing) a re-brightening phase in 2012 we extend the X-ray 
luminosity range. Our 2012 data show that the source ap¬ 
pears to lie on same steep track and at 5xl0 34 ergs~ 1 the 
source is still radio under-luminous compared to other X- 
ray binaries at comparable Lx (Figure (9]). The relation we 
observe also holds over ~2.5 orders of magnitude in radio 
luminosity and does not appear to be random scatter about 
the standard correlation. Gi ven the low inclinat ion of this 
system (between 4° and 15°: lRussell et ahll2014bl) it is plau¬ 
sible that variable Doppler boosting affects the correlation. 
We now investigate this possibility. 

Assuming random inclina t ion a ngles and highly vari¬ 
able boosting, ISoIeri fc Fende/i ( 2011 ) tested a model to de¬ 
termine if variable boosting could account for the scatter of 
the radio/X-ray correlation for the available sample of black 


© 2015 RAS, MNRAS OOO.ITIIHI 



















































Radio monitoring of the hard state jets in MAXI J1836—19f 13 


holes. With their model, they were able to qualitatively re¬ 
produce both the scatter of the correlation, as well as the 
‘radio-quiet’ lower branch but noted that the model had its 
limitations. Here, we use a similar approach to test if the 
steep radio/X-ray correlation could be explained by vari¬ 
able relativistic beaming. While recent work has cast doubt 
over the universality of the two well-defin ed tracks for the 
radio/X-ray correlation dGallo et al.ll2014h . we assume that 
the intrinsic radio luminosity of MAXI J1836—194 lies on the 
lower canonical hard-state tracks. We then determine what 
boosting would be required to account for the observed ra¬ 
dio luminosity, testing whether a variable T can plausibly 
account for the steeper than usual correlation. For the sys¬ 
tem to lie on the upper track, the radio emission at low 
X-ray luminosities would need to be significantly deboosted 
(while being boosted at higher X-ray luminosities; Figure[9]). 
Deboosted radio e mission has been observed in LMXBs (e.g. 
XTE J1752—223; I Yang et al .11201(1 ) but is unlikely to have 
affected our observations due to the low inclination angle of 
the system. 

We assume that the X-ray emission remains unbeamed 
(but see lMarkoff et af]|2005l : [Russell et al.ll2oiol l and relate 
the observed beamed radio luminosity, Lr, 0 , to the intrinsic 
radio luminosity, Lr,:, by 

fR,o = W"-“, (5) 

where 5 is the Doppl er fac tor, n=2 (for a continuously re¬ 
plenished jet; lFendeij|2006l l. and a is the spectral index of 
the radio emission. The Doppler factor is a function of F 
where <5 = r _1 (l — /3cos#) -1 , /3 is the velocity of the jet 
emission as a fraction of the speed of light, and 9 is the in¬ 
clination of the system; because r = (1 — (I 2 ) -1 / 2 , we are 
then able to calculate the Lorentz factor required to shift 
the modelled intrinsic radio luminosity to tha t observed. 

_ If we assume an inclination angle of 10° llRussell et al.1 

l2014blf . a source distance of 10 kpc (iRussell et al.ll2014al l 
and a jet axis perpendicular to the plane of the binary, then 
for Ln,i to lie on the lower track of the radio/X-ray correla¬ 
tion T must change by a factor of ~3-4 during the outburst 
(from F «1 at the lowest X-ray luminosity to T ~3-4 at the 
highest; Figure HOT) . This result is sensitive to the normalisa¬ 
tion of the assumed intrinsic radio luminosity, therefore we 
have tested a range of normalisations that account for the 
observed scatter about the lower track. The maximum values 
of F are higher th an previously thou ght for both a compact 
jet (where F +2; iFender et al]|2004l f a nd the transient je t 
of GRS 1915+105 (where 1.6+ T +1.9: [Reid et al.ll2014al l. 
However, due to the low inclination of MAXI J1836—194, 
the values we estimate are not so high that beaming can be 
ruled out as an explanation for the steep radio/X-ray cor¬ 
relation. Fitting r against X-ray luminosity (as a fraction 
of the Eddington luminosity for a 10 Mq black hole) gives 
log(F) oc reZ x , where l x = L x /L ecM and re « 100-120. With 
this result, the bulk Lorentz factor become s larger than ~1 
above 0.05 percent LecM- For comparison, ISoleri fc Fender! 
(1201 lT) found F becomes > 1 above 0.1 percent Z/Edd- This 
level of boosting would be expected to have consequences on 
the observed radio spectrum, which coupled with the radio 
spectral index may then provide some insight into physical 
parameters of the jet (such as the opening angle). However, 
a full modelling analysis would be required to determine the 
jet parameters which is beyond the scope of this paper. 


At a source distance of <8 kpc, the radio emission at 
low X-ray luminosities falls well below the canonical rela¬ 
tion (Figure [21) and could not be explained unless it was 
deboosted. Due to the low inclination angle of the system, 
we do not expect deboosting. Therefore, if we assume that 
the steep radio/X-ray correlation is due to a variable Lorentz 
factor, MAXI J1836—194 must be at a distance greater than 
8 kpc. 

Our results show that relativistic beaming can 
account for the observed radio/X-ray correlation in 
MAXI J1836—194. However, similar to the results of 
ISoleri fc Fender! d201lf) , if we applied our boosting argument 
to the full sample of LMXBs we would see the correlation 
break down at X-ray luminosities > 0.5 percent Eddington. 
Following an outburst, black hole LMXBs generally tran¬ 
sition back to the hard state at luminositie s between 0.3 
and 3 percent Eddin gton dKale mci et al.ll2013h , with a mean 
value of ~2 percent dMaccaron d M), demonstrating that 
X-ray luminosities of > 0.5 p ercent Eddington are ex pected 
during the hard state decay. iHeinz fc Merloni d2004l l anal¬ 
ysed the radio/X-ray relation of GX 339—4 and V404 Cyg 
and argued that similar Lorentz factors are expected be¬ 
tween sources. We therefore conclude that while a variable 
Lorentz factor can explain the unusually steep radio/X-ray 
correlation in MAXI J1836—194, it cannot be universally ap¬ 
plied to the broad sample of black hole binaries. 

It is possible that the steep corr elation is a temporar y 
deviation from the standard track dCorbel et al.l l2013blb 
However, when we include data from the 2012 re-flare, we see 
the source appears to lie on the same steep radio/X-ray track 
and the steep correlation and does not immediately rejoin 
the standard track. This system could simply produce an in- 
trinsically steep c orrelation (~2cr fro m the steeper trac k of 
ICoriat et al.ll2()ill . or 4a from that of Gallo et al.ll201 j . but 
see recent results from lGallo et al.ll2014 l. however, this index 
cannot be explained by models of X-ray emission originating 
as inverse Compton from a corona, or being dominated by 
optically-thin synchrotron emission from the compact jet. 


5 CONCLUSIONS 

We have presented the full VLA, ATCA and VLBA moni¬ 
toring of MAXI J1836—194 during its 2011 outburst. With 
our intensive radio coverage of the system we have observed 
the evolution of the polarised compact jet. We find that the 
compact jet was linearly polarised at a few percent during 
the early, bright phase of the outburst. At later times we 
were unable to place any constraint on the linear polarisa¬ 
tion from the source. We determined the VLBI jet axis and 
found an alignment with the EVPA (which was calculated 
from the measured Stokes parameters). This implies that the 
magnetic field was perpendicular to the compact jet axis, 
possibly due to a helical structure from the rotation of the 
accretion disk or compression due to internal shocks within 
the compact jet. If this alignment is common for all LMXBs 
then high resolution VLBI observations may not be required 
to determine the orientation of the jet. However, further 
studies are required to gauge whether the EVPA is always 
aligned with the jet axis. Calculating the peculiar velocity 
from our astrometric observations, we find that the black 
hole likely required an asymmetric natal kick during for- 
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mation to account for its high peculiar velocity. During the 
decay phase of the outburst we observed an unusually steep 
radio/X-ray correlation where Lr oc &/ 8±0 - 2 5 much steeper 
than other black hole LMXB systems. We test whether the 
steep correlation may be due to a variable jet Lorentz fac¬ 
tor. If we assume that the correlation agrees with the radio 
quiet radio/X-ray correlation track, the Lorentz factor needs 
to vary from T ~1 at low X-ray luminosities to L ~3-4 at 
high X-ray luminosities to produce the observed correlation. 
However, this relation is unlikely to explain to the full sam¬ 
ple of black hole LMXBs, and it is possible that we are ob¬ 
serving a temporary deviation to the standard correlation. 
Interestingly, data from a re-flare event in 2012 show that the 
source appears to once again lie on the unusually steep track 
followed during the 2011 outburst. If the steep radio/X-ray 
correlation is due only to a variable Lorentz factor, our data 
require that the source lies at a distance greater than 8kpc. 
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Table 1: continued VLA flux densities of MAXI J1836—194. a is the radio 
spectral index. Quoted 1 <t errors are uncertainties to the fitted source 
model. 


Date 

(UT) 

MJD 

Frequency 

(GHz) 

Flux density 
(mJy) 

Q 

(mJy beam x ) 

U 

(mJy beam -1 ) 

a 

2012 Sep 03 

55807.1 

4.60 

27.2±0.3 

0.24±0.05 

-0.29±0.06 

0.64A0.03 



7.90 

38.4A0.5 

0.67±0.06 

-0.49±0.07 


2011 Sep 05 

55809.1 

5.00 

23.2±0.3 

0.08±0.02 

-0.08±0.02 

0.84A0.04 



7.45 

32.5±0.4 

0.18±0.02 

-0.16±0.02 


2011 Sep 07 

55811.2 

1.50 

18.1±0.4 

- 

- 

0.51i0.01 



5.00 

28.6±0.7 

0.25T0.07 

-0.13±0.07 




7.45 

39.3±0.9 

0.58±0.09 

-0.08±0.08 




20.80 

68±2 

- 

- 




25.90 

79±2 

- 

- 




32.02 

88±5 

- 

- 




41.00 

87±8 

- 

- 


2011 Sep 09 

55813.2 

1.60 

24.9±0.5 

- 

- 

0.43A0.01 



5.00 

33.7±0.6 

-0.00±0.06 

-0.15±0.05 




7.45 

41.4±0.9 

0.29±0.05 

-0.49±0.05 




20.80 

66±2 

- 

- 




25.90 

73±2 

- 

- 




32.02 

95±5 

- 

- 




41.00 

108±6 

- 

- 


2011 Sep 12 

55816.1 

1.60 

36.7±0.6 

- 

- 

0.10±0.01 



5.00 

36.6±0.8 

0.93±0.05 

-0.94±0.05 




7.45 

39.5±0.8 

0.72±0.05 

-0.74±0.04 




20.80 

47±2 

- 

- 




25.90 

48±2 

- 

- 




32.02 

51±3 

- 

- 




41.00 

55±3 

- 

- 


2011 Sep 13 

55817.0 

1.60 

29.0±0.5 

- 

- 

0.14i0.01 



5.00 

31.0±0.4 

0.72±0.06 

-0.92±0.06 




7.45 

33.1±0.4 

0.62±0.05 

-0.86±0.05 




20.80 

40.7±1 

- 

- 




25.90 

44±2 

- 

- 




32.02 

46±3 

- 

- 




41.00 

46A3 

- 

- 


2011 Sep 18 

55822.0 

1.50 

37il 

- 

- 

0.13i0.01 



5.00 

34.7±0.4 

0.85±0.06 

-0.85±0.06 




7.45 

36.9±0.9 

0.62±0.06 

-0.59±0.06 




20.80 

44±1 

- 

- 




25.90 

46±2 

- 

- 




32.02 

50±2 

- 

- 




41.00 

441±3 

- 

- 


2011 Sep 18 

55822.9 

1.56 

33A2 

- 

- 

0.14A0.01 



5.00 

29.7±0.4 

0.72±0.07 

-0.75±0.07 




7.45 

31.3±0.4 

0.58±0.06 

-0.56i0.06 




20.80 

35.5±0.7 

- 

- 




25.90 

38±2 

- 

- 




32.02 

44±2 

- 

- 




41.00 

37±2 

- 

- 


2011 Sep 23 

55827.061 

5.00 

32.7±0.8 

- 

- 

0.34±0.02 



7.45 

34.0±0.8 

- 

- 




20.80 

51±2 

- 

- 




25.90 

54±2 

- 

- 




32.02 

62±3 

- 

- 




41.00 

66A4 

- 

- 


2011 Sep 26 

55830.9 

5.00 

14.1±0.4 

- 

- 

0.51A0.02 



7.45 

14.8±0.4 

- 

- 


2011 Sep 26 

55830.9 

20.80 

27±1 

- 

- 




25.90 

30il 

- 

- 

0.51i0.02 
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Date 

(UT) 

MJD 

Frequency 

(GHz) 

Flux density 

(mJy) 

Q 

(mJy beam 

U 

(mJy beam -1 ) 

a 



31.52 

36±2 

- 

- 




37.52 

40±2 

- 

- 


2011 Sep 30 

55834.0 

5.26 

11.3±0.3 

- 

- 

0.57±0.02 



7.45 

11.7±0.3 

- 

- 




20.80 

23.5±0.9 

- 

- 




25.90 

26±1 

- 

- 




31.52 

30±2 

- 

- 




37.52 

32±2 

- 

- 


2011 Oct 07 

55841.9 

5.26 

8±0.3 

- 

- 

0.59±0.03 



7.45 

8±0.3 

- 

- 




20.80 

17.5±0.9 

- 

- 




25.90 

19.3±0.9 

- 

- 




31.52 

22±2 

- 

- 




37.52 

25±2 

- 

- 


2011 Oct 12 

55846.0 

5.26 

5.0±0.3 

- 

- 

0.55±0.03 



7.45 

5.7±0.2 

- 

- 




20.80 

9.8±0.5 

- 

- 




25.90 

12.2±0.7 

- 

- 




31.52 

Mil 

- 

- 




37.52 

14±1 

- 

- 


2011 Oct 22 

55856.9 

5.26 

3.8±0.2 

- 

- 

0.25±0.04 



7.45 

4.0±0.2 

- 

- 




20.80 

5.3±0.4 

- 

- 




25.90 

5.6±0.4 

- 

- 




31.52 

5.8±0.6 

- 

- 




37.52 

5.7±0.8 

- 

- 


2011 Nov 01 

55866.9 

5.26 

1.5±0.1 

- 

- 

0.33±0.06 



7.45 

1.46±0.08 

- 

- 




20.80 

1.9±0.2 

- 

- 




25.90 

2.6±0.4 

- 

- 




31.52 

2.6±0.4 

- 

- 




37.52 

3.2±0.7 

- 

- 


2011 Nov 11 

55876.8 

5.26 

0.8±0.2 

- 

- 

0.4±0.1 



7.45 

0.8±0.1 

- 

- 




20.80 

1.0±0.1 

- 

- 




25.90 

1.5±0.2 

- 

- 


2011 Nov 18 

55883.9 

5.26 

0.33±0.03 

- 

- 

0.5±0.1 



7.45 

0.31±0.06 

- 

- 




20.80 

0.6±0.1 

- 

- 




25.90 

0.8±0.1 

- 

- 


2011 Dec 03 

55898.8 

5.26 

0.19±0.04 

- 

- 

0.7±0.1 



7.45 

0.25±0.02 

- 

- 




20.80 

0.49±0.07 

- 

- 




25.90 

0.5±0.1 

- 

- 
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